(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.4' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     78478,       2643]
NotebookOptionsPosition[     73094,       2487]
NotebookOutlinePosition[     75663,       2549]
CellTagsIndexPosition[     75577,       2544]
WindowTitle->DoPolarizationSums
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[GridBox[{
   {Cell["FEYN CALC SYMBOL", "PacletNameCell"], Cell[TextData[{
     Cell[BoxData[
      ActionMenuBox[
       FrameBox[
        InterpretationBox[Cell[TextData[{
          "See Also",
          StyleBox[" \[FilledDownTriangle]", "AnchorBarArrow",
           StripOnInput->False]
         }]],
         TextCell[
          Row[{"See Also", 
            Style[" \[FilledDownTriangle]", "AnchorBarArrow"]}]]],
        StripOnInput->False], {"\<\"Polariazation\"\>":>
       Documentation`HelpLookup[
       "paclet:FeynCalc/ref/Polariazation"], "\<\"Uncontract\"\>":>
       Documentation`HelpLookup["paclet:FeynCalc/ref/Uncontract"]},
       Appearance->None,
       MenuAppearance->Automatic,
       MenuStyle->"SeeAlso"]],
      LineSpacing->{1.4, 0}],
     "\[ThickSpace]\[ThickSpace]\[ThickSpace]\[ThickSpace]\[ThickSpace]\
\[ThickSpace]",
     Cell[BoxData[
      ActionMenuBox[
       FrameBox[
        InterpretationBox[Cell[TextData[{
          "URL",
          StyleBox[" \[FilledDownTriangle]", "AnchorBarArrow",
           StripOnInput->False]
         }]],
         TextCell[
          Row[{"URL", 
            Style[" \[FilledDownTriangle]", "AnchorBarArrow"]}]]],
        StripOnInput->False], {"\<\"FeynCalc/ref/DoPolarizationSums\"\>":>
       None, "\<\"Copy Wolfram Documentation Center URL\"\>":>
       Module[{DocumentationSearch`Private`nb$}, 
        DocumentationSearch`Private`nb$ = NotebookPut[
           Notebook[{
             Cell["FeynCalc/ref/DoPolarizationSums"]}, Visible -> False]]; 
        SelectionMove[DocumentationSearch`Private`nb$, All, Notebook]; 
        FrontEndTokenExecute[DocumentationSearch`Private`nb$, "Copy"]; 
        NotebookClose[DocumentationSearch`Private`nb$]; Null], 
       Delimiter, "\<\"Copy web URL\"\>":>
       Module[{DocumentationSearch`Private`nb$}, 
        DocumentationSearch`Private`nb$ = NotebookPut[
           Notebook[{
             Cell[
              BoxData[
               MakeBoxes[
                Hyperlink[
                "http://reference.wolfram.com/language/FeynCalc/ref/\
DoPolarizationSums.html"], StandardForm]], "Input", TextClipboardType -> 
              "PlainText"]}, Visible -> False]]; 
        SelectionMove[DocumentationSearch`Private`nb$, All, Notebook]; 
        FrontEndTokenExecute[DocumentationSearch`Private`nb$, "Copy"]; 
        NotebookClose[DocumentationSearch`Private`nb$]; 
        Null], "\<\"Go to web URL\"\>":>FrontEndExecute[{
         NotebookLocate[{
           URL[
            StringJoin[
             If[
              TrueQ[DocumentationBuild`Make`Private`wsmlinkQ$73676], 
              "http://reference.wolfram.com/system-modeler/", 
              "http://reference.wolfram.com/language/"], 
             "FeynCalc/ref/DoPolarizationSums", ".html"]], None}]}]},
       Appearance->None,
       MenuAppearance->Automatic,
       MenuStyle->"URLMenu"]],
      LineSpacing->{1.4, 0}]
    }], "AnchorBar"]}
  }]], "AnchorBarGrid",
 GridBoxOptions->{GridBoxItemSize->{"Columns" -> {
     Scaled[0.65], {
      Scaled[0.34]}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
   "RowsIndexed" -> {}}},
 CellID->1],

Cell[TextData[{
 Cell["DoPolarizationSums", "ObjectName"],
 Cell[BoxData[
  InterpretationBox[
   StyleBox[
    GraphicsBox[{},
     BaselinePosition->Baseline,
     ImageSize->{8, 0}],
    CacheGraphics->False],
   Spacer[8]]]],
 Cell[BoxData[""], "ObjectNameTranslation"]
}], "ObjectNameGrid"],

Cell[CellGroupData[{

Cell[BoxData[GridBox[{
   {"", Cell[TextData[{
     Cell[BoxData[
      RowBox[{"DoPolarizationSums", "[", 
       RowBox[{"exp", ",", "k", ",", "n"}], "]"}]], "InlineFormula"],
     " \[LineSeparator]sums over physical (transverse) polarizations of an \
external massless vector boson with momentum ",
     Cell[BoxData["k"], "InlineFormula"],
     ", where ",
     Cell[BoxData["n"], "InlineFormula"],
     " is an auxiliary 4-vector from the gauge-dependent polarization sum \
formula."
    }]]},
   {"", Cell[TextData[{
     Cell[BoxData[
      RowBox[{"DoPolarizationSums", "[", 
       RowBox[{"exp", ",", "k", ",", "0"}], "]"}]], "InlineFormula"],
     " \[LineSeparator]replaces the polarization sum of an external massless \
vector boson with momentum ",
     Cell[BoxData["k"], "InlineFormula"],
     " by ",
     Cell[BoxData[
      RowBox[{"-", 
       SuperscriptBox["g", "\[Mu]\[Nu]"]}]], "InlineFormula"],
     ". This corresponds to the summation over all four polarizations, \
including the unphysical ones."
    }]]},
   {"", Cell[TextData[{
     Cell[BoxData[
      RowBox[{"DoPolarizationSums", "[", 
       RowBox[{"exp", ",", "k"}], "]"}]], "InlineFormula"],
     " \[LineSeparator]sums over the three polarizations of an external \
massive vector boson with momentum ",
     Cell[BoxData["k"], "InlineFormula"],
     " and mass ",
     Cell[BoxData[
      SuperscriptBox["k", "2"]], "InlineFormula"],
     "."
    }]]}
  }]], "Usage",
 GridBoxOptions->{
 GridBoxBackground->{
  "Columns" -> {{None}}, "ColumnsIndexed" -> {}, 
   "Rows" -> {None, None, {None}}, "RowsIndexed" -> {}}},
 CellID->2010081510],

Cell[CellGroupData[{

Cell[TextData[Cell[BoxData[
 ButtonBox[Cell[TextData[{
   Cell[BoxData[
    InterpretationBox[
     StyleBox[
      GraphicsBox[{},
       BaselinePosition->Baseline,
       ImageSize->{6, 0}],
      CacheGraphics->False],
     Spacer[6]]]],
   "Details"
  }], "NotesFrameText"],
  Appearance->{Automatic, None},
  BaseStyle->None,
  ButtonFunction:>(FrontEndExecute[{
     FrontEnd`SelectionMove[
      FrontEnd`SelectedNotebook[], All, ButtonCell], 
     FrontEndToken["OpenCloseGroup"], 
     FrontEnd`SelectionMove[
      FrontEnd`SelectedNotebook[], After, CellContents]}]& ),
  Evaluator->None,
  Method->"Preemptive"]]]], "NotesSection",
 WholeCellGroupOpener->True,
 CellGroupingRules->{"SectionGrouping", 50},
 CellID->515142955],

Cell["", "SectionHeaderSpacer"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Options", "[", "DoPolarizationSums", "]"}]], "Input",
 CellLabel->"In[39]:="],

Cell[BoxData[
 FormBox[
  RowBox[{"{", 
   RowBox[{
    RowBox[{"Contract", "\[Rule]", "True"}], ",", 
    RowBox[{"ExtraFactor", "\[Rule]", "1"}], ",", 
    RowBox[{"FeynCalcExternal", "\[Rule]", "False"}], ",", 
    RowBox[{"FeynCalcInternal", "\[Rule]", "False"}], ",", 
    RowBox[{"GaugeTrickN", "\[Rule]", "2"}], ",", 
    RowBox[{"Head", "\[Rule]", "Identity"}], ",", 
    RowBox[{"VirtualBoson", "\[Rule]", "False"}]}], "}"}], 
  TraditionalForm]], "Output",
 ImageSize->{518, 35},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[39]="]
}, Open  ]]
}, Closed]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 Cell[BoxData[
  InterpretationBox[
   StyleBox[
    GraphicsBox[{},
     BaselinePosition->Baseline,
     ImageSize->{6, 0}],
    CacheGraphics->False],
   Spacer[6]]]],
 "Examples",
 "\[NonBreakingSpace]\[NonBreakingSpace]",
 Cell["(5)", "ExampleCount"]
}], "PrimaryExamplesSection",
 WholeCellGroupOpener->True,
 CellTags->"PrimaryExamplesSection",
 CellID->95617754],

Cell[CellGroupData[{

Cell[TextData[{
 "Basic Examples",
 "\[NonBreakingSpace]\[NonBreakingSpace]",
 Cell["(5)", "ExampleCount"]
}], "ExampleSection", "ExampleSection",
 WholeCellGroupOpener->True,
 CellID->1253010737],

Cell["\<\
The standard formula for massless vector bosons is valid for all types of the \
corresponding particles, including gluons.\
\>", "Notes"],

Cell[BoxData[{
 RowBox[{"FCClearScalarProducts", "[", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "p", "]"}], "=", "0"}], ";"}]}], "Input",
 CellLabel->"In[1]:="],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Pair", "[", 
   RowBox[{
    RowBox[{"LorentzIndex", "[", "\[Mu]", "]"}], ",", 
    RowBox[{"Momentum", "[", 
     RowBox[{"Polarization", "[", 
      RowBox[{"p", ",", 
       RowBox[{"-", "\[ImaginaryI]"}]}], "]"}], "]"}]}], "]"}], " ", 
  RowBox[{"Pair", "[", 
   RowBox[{
    RowBox[{"LorentzIndex", "[", "\[Nu]", "]"}], ",", 
    RowBox[{"Momentum", "[", 
     RowBox[{"Polarization", "[", 
      RowBox[{"p", ",", "\[ImaginaryI]"}], "]"}], "]"}]}], "]"}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{
    SuperscriptBox[
     SuperscriptBox[
      OverscriptBox["\[CurlyEpsilon]", "_"], "*"], 
     FormBox[
      FormBox["\[Mu]",
       TraditionalForm],
      TraditionalForm]], "(", 
    FormBox["p",
     TraditionalForm], ")"}], " ", 
   RowBox[{
    SuperscriptBox[
     OverscriptBox["\[CurlyEpsilon]", "_"], 
     FormBox[
      FormBox["\[Nu]",
       TraditionalForm],
      TraditionalForm]], "(", 
    FormBox["p",
     TraditionalForm], ")"}]}], TraditionalForm]], "Output",
 ImageSize->{101, 16},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[2]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DoPolarizationSums", "[", 
  RowBox[{"%", ",", "p", ",", "n"}], "]"}]], "Input",
 CellLabel->"In[3]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      SuperscriptBox[
       FormBox[
        OverscriptBox[
         FormBox["n",
          TraditionalForm], "_"],
        TraditionalForm], "2"], " ", 
      SuperscriptBox[
       FormBox[
        OverscriptBox[
         FormBox["p",
          TraditionalForm], "_"],
        TraditionalForm], 
       FormBox[
        FormBox["\[Mu]",
         TraditionalForm],
        TraditionalForm]], " ", 
      SuperscriptBox[
       FormBox[
        OverscriptBox[
         FormBox["p",
          TraditionalForm], "_"],
        TraditionalForm], 
       FormBox[
        FormBox["\[Nu]",
         TraditionalForm],
        TraditionalForm]]}], 
     SuperscriptBox[
      RowBox[{
       FormBox["\<\"(\"\>",
        TraditionalForm], 
       FormBox[
        RowBox[{
         FormBox[
          OverscriptBox[
           FormBox["n",
            TraditionalForm], "_"],
          TraditionalForm], 
         FormBox["\<\"\[CenterDot]\"\>",
          TraditionalForm], 
         FormBox[
          OverscriptBox[
           FormBox["p",
            TraditionalForm], "_"],
          TraditionalForm]}],
        TraditionalForm], 
       FormBox["\<\")\"\>",
        TraditionalForm]}], "2"]]}], "-", 
   SuperscriptBox[
    OverscriptBox["g", "_"], 
    RowBox[{
     FormBox[
      FormBox["\[Mu]",
       TraditionalForm],
      TraditionalForm], 
     FormBox[
      FormBox["\[Nu]",
       TraditionalForm],
      TraditionalForm]}]], "+", 
   FractionBox[
    RowBox[{
     SuperscriptBox[
      FormBox[
       OverscriptBox[
        FormBox["n",
         TraditionalForm], "_"],
       TraditionalForm], 
      FormBox[
       FormBox["\[Nu]",
        TraditionalForm],
       TraditionalForm]], " ", 
     SuperscriptBox[
      FormBox[
       OverscriptBox[
        FormBox["p",
         TraditionalForm], "_"],
       TraditionalForm], 
      FormBox[
       FormBox["\[Mu]",
        TraditionalForm],
       TraditionalForm]]}], 
    RowBox[{
     FormBox[
      OverscriptBox[
       FormBox["n",
        TraditionalForm], "_"],
      TraditionalForm], 
     FormBox["\<\"\[CenterDot]\"\>",
      TraditionalForm], 
     FormBox[
      OverscriptBox[
       FormBox["p",
        TraditionalForm], "_"],
      TraditionalForm]}]], "+", 
   FractionBox[
    RowBox[{
     SuperscriptBox[
      FormBox[
       OverscriptBox[
        FormBox["n",
         TraditionalForm], "_"],
       TraditionalForm], 
      FormBox[
       FormBox["\[Mu]",
        TraditionalForm],
       TraditionalForm]], " ", 
     SuperscriptBox[
      FormBox[
       OverscriptBox[
        FormBox["p",
         TraditionalForm], "_"],
       TraditionalForm], 
      FormBox[
       FormBox["\[Nu]",
        TraditionalForm],
       TraditionalForm]]}], 
    RowBox[{
     FormBox[
      OverscriptBox[
       FormBox["n",
        TraditionalForm], "_"],
      TraditionalForm], 
     FormBox["\<\"\[CenterDot]\"\>",
      TraditionalForm], 
     FormBox[
      OverscriptBox[
       FormBox["p",
        TraditionalForm], "_"],
      TraditionalForm]}]]}], TraditionalForm]], "Output",
 ImageSize->{277, 43},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[3]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[Cell["\t", "ExampleDelimiter"],
  $Line = 0; Null]], "ExampleDelimiter"],

Cell["\<\
In QED the gauge invariance ensures the cancellation of unphysical \
polarizations so that for photons one can also employ the simpler replacement \
with the metric tensor.\
\>", "Notes"],

Cell[BoxData[{
 RowBox[{"FCClearScalarProducts", "[", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "p", "]"}], "=", "0"}], ";"}]}], "Input",
 CellLabel->"In[1]:="],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Pair", "[", 
   RowBox[{
    RowBox[{"LorentzIndex", "[", "\[Mu]", "]"}], ",", 
    RowBox[{"Momentum", "[", 
     RowBox[{"Polarization", "[", 
      RowBox[{"p", ",", 
       RowBox[{"-", "\[ImaginaryI]"}]}], "]"}], "]"}]}], "]"}], " ", 
  RowBox[{"Pair", "[", 
   RowBox[{
    RowBox[{"LorentzIndex", "[", "\[Nu]", "]"}], ",", 
    RowBox[{"Momentum", "[", 
     RowBox[{"Polarization", "[", 
      RowBox[{"p", ",", "\[ImaginaryI]"}], "]"}], "]"}]}], "]"}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{
    SuperscriptBox[
     SuperscriptBox[
      OverscriptBox["\[CurlyEpsilon]", "_"], "*"], 
     FormBox[
      FormBox["\[Mu]",
       TraditionalForm],
      TraditionalForm]], "(", 
    FormBox["p",
     TraditionalForm], ")"}], " ", 
   RowBox[{
    SuperscriptBox[
     OverscriptBox["\[CurlyEpsilon]", "_"], 
     FormBox[
      FormBox["\[Nu]",
       TraditionalForm],
      TraditionalForm]], "(", 
    FormBox["p",
     TraditionalForm], ")"}]}], TraditionalForm]], "Output",
 ImageSize->{101, 16},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[2]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DoPolarizationSums", "[", 
  RowBox[{"%", ",", "p", ",", "0"}], "]"}]], "Input",
 CellLabel->"In[3]:="],

Cell[BoxData[
 FormBox[
  RowBox[{"-", 
   SuperscriptBox[
    OverscriptBox["g", "_"], 
    RowBox[{
     FormBox[
      FormBox["\[Mu]",
       TraditionalForm],
      TraditionalForm], 
     FormBox[
      FormBox["\[Nu]",
       TraditionalForm],
      TraditionalForm]}]]}], TraditionalForm]], "Output",
 ImageSize->{48, 16},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[3]="]
}, Open  ]],

Cell["\<\
You can also use this trick in QCD, provided that unphysical degress of \
freedom are subtracted using ghosts at a later stage. Notice that in this \
case you should not make the polarization vectors transverse using the \
Transversality option. Furthermore, the averaging over the polarizations of \
the initial gluons must done on the physical amplitude squared, i.e. after \
the ghost contributions have been subtracted.\
\>", "Notes"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[Cell["\t", "ExampleDelimiter"],
  $Line = 0; Null]], "ExampleDelimiter"],

Cell["\<\
Massive vector bosons (e.g. W or Z) have 3 degrees of freedom and require no \
auxiliary vector.\
\>", "Notes"],

Cell[BoxData[{
 RowBox[{"FCClearScalarProducts", "[", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "p", "]"}], "=", 
   RowBox[{"m", "^", "2"}]}], ";"}]}], "Input",
 CellLabel->"In[1]:="],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Pair", "[", 
   RowBox[{
    RowBox[{"LorentzIndex", "[", "\[Mu]", "]"}], ",", 
    RowBox[{"Momentum", "[", 
     RowBox[{"Polarization", "[", 
      RowBox[{"p", ",", 
       RowBox[{"-", "\[ImaginaryI]"}]}], "]"}], "]"}]}], "]"}], " ", 
  RowBox[{"Pair", "[", 
   RowBox[{
    RowBox[{"LorentzIndex", "[", "\[Nu]", "]"}], ",", 
    RowBox[{"Momentum", "[", 
     RowBox[{"Polarization", "[", 
      RowBox[{"p", ",", "\[ImaginaryI]"}], "]"}], "]"}]}], "]"}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{
    SuperscriptBox[
     SuperscriptBox[
      OverscriptBox["\[CurlyEpsilon]", "_"], "*"], 
     FormBox[
      FormBox["\[Mu]",
       TraditionalForm],
      TraditionalForm]], "(", 
    FormBox["p",
     TraditionalForm], ")"}], " ", 
   RowBox[{
    SuperscriptBox[
     OverscriptBox["\[CurlyEpsilon]", "_"], 
     FormBox[
      FormBox["\[Nu]",
       TraditionalForm],
      TraditionalForm]], "(", 
    FormBox["p",
     TraditionalForm], ")"}]}], TraditionalForm]], "Output",
 ImageSize->{101, 16},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[2]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DoPolarizationSums", "[", 
  RowBox[{"%", ",", "p"}], "]"}]], "Input",
 CellLabel->"In[3]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   FractionBox[
    RowBox[{
     SuperscriptBox[
      FormBox[
       OverscriptBox[
        FormBox["p",
         TraditionalForm], "_"],
       TraditionalForm], 
      FormBox[
       FormBox["\[Mu]",
        TraditionalForm],
       TraditionalForm]], " ", 
     SuperscriptBox[
      FormBox[
       OverscriptBox[
        FormBox["p",
         TraditionalForm], "_"],
       TraditionalForm], 
      FormBox[
       FormBox["\[Nu]",
        TraditionalForm],
       TraditionalForm]]}], 
    SuperscriptBox["m", "2"]], "-", 
   SuperscriptBox[
    OverscriptBox["g", "_"], 
    RowBox[{
     FormBox[
      FormBox["\[Mu]",
       TraditionalForm],
      TraditionalForm], 
     FormBox[
      FormBox["\[Nu]",
       TraditionalForm],
      TraditionalForm]}]]}], TraditionalForm]], "Output",
 ImageSize->{106, 39},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[3]="]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[Cell["\t", "ExampleDelimiter"],
  $Line = 0; Null]], "ExampleDelimiter"],

Cell[TextData[{
 "A more realistic example of summing over the polarizations of the photons \
in ",
 Cell[BoxData[
  RowBox[{
   RowBox[{
    SuperscriptBox["e", "+"], 
    SuperscriptBox["e", "-"]}], "\[Rule]", " ", "\[Gamma]\[Gamma]"}]], 
  "InlineFormula"]
}], "Notes"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ClearAll", "[", 
   RowBox[{"s", ",", "t", ",", "u"}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"FCClearScalarProducts", "[", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "k1", "]"}], "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "k2", "]"}], "=", "0"}], ";"}]}], "Input",
 CellLabel->"In[1]:="],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"Spinor", "[", 
           RowBox[{
            RowBox[{"Momentum", "[", "p1", "]"}], ",", "0", ",", "1"}], "]"}],
           ".", 
          RowBox[{"GS", "[", 
           RowBox[{"Polarization", "[", 
            RowBox[{"k1", ",", "I", ",", 
             RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], 
          ".", 
          RowBox[{"GS", "[", 
           RowBox[{"k2", "-", "p2"}], "]"}], ".", 
          RowBox[{"GS", "[", 
           RowBox[{"Polarization", "[", 
            RowBox[{"k2", ",", "I", ",", 
             RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], 
          ".", 
          RowBox[{"Spinor", "[", 
           RowBox[{
            RowBox[{"-", 
             RowBox[{"Momentum", "[", "p2", "]"}]}], ",", "0", ",", "1"}], 
           "]"}]}], "*", 
         RowBox[{
          RowBox[{"SMP", "[", "\"\<e\>\"", "]"}], "^", "2"}]}], ")"}], "/", 
       "t"}], ")"}]}], "-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"Spinor", "[", 
         RowBox[{
          RowBox[{"Momentum", "[", "p1", "]"}], ",", "0", ",", "1"}], "]"}], 
        ".", 
        RowBox[{"GS", "[", 
         RowBox[{"Polarization", "[", 
          RowBox[{"k2", ",", "I", ",", 
           RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], ".", 
        RowBox[{"GS", "[", 
         RowBox[{"k1", "-", "p2"}], "]"}], ".", 
        RowBox[{"GS", "[", 
         RowBox[{"Polarization", "[", 
          RowBox[{"k1", ",", "I", ",", 
           RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], ".", 
        RowBox[{"Spinor", "[", 
         RowBox[{
          RowBox[{"-", 
           RowBox[{"Momentum", "[", "p2", "]"}]}], ",", "0", ",", "1"}], 
         "]"}]}], "*", 
       RowBox[{
        RowBox[{"SMP", "[", "\"\<e\>\"", "]"}], "^", "2"}]}], ")"}], "/", 
     "u"}]}], ")"}], "*", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"Spinor", "[", 
           RowBox[{
            RowBox[{"-", 
             RowBox[{"Momentum", "[", "p2", "]"}]}], ",", "0", ",", "1"}], 
           "]"}], ".", 
          RowBox[{"GS", "[", 
           RowBox[{"Polarization", "[", 
            RowBox[{"k1", ",", 
             RowBox[{"-", "I"}], ",", 
             RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], 
          ".", 
          RowBox[{"GS", "[", 
           RowBox[{"k1", "-", "p2"}], "]"}], ".", 
          RowBox[{"GS", "[", 
           RowBox[{"Polarization", "[", 
            RowBox[{"k2", ",", 
             RowBox[{"-", "I"}], ",", 
             RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], 
          ".", 
          RowBox[{"Spinor", "[", 
           RowBox[{
            RowBox[{"Momentum", "[", "p1", "]"}], ",", "0", ",", "1"}], 
           "]"}]}], "*", 
         RowBox[{
          RowBox[{"SMP", "[", "\"\<e\>\"", "]"}], "^", "2"}]}], ")"}], "/", 
       "u"}], ")"}]}], "-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"Spinor", "[", 
         RowBox[{
          RowBox[{"-", 
           RowBox[{"Momentum", "[", "p2", "]"}]}], ",", "0", ",", "1"}], 
         "]"}], ".", 
        RowBox[{"GS", "[", 
         RowBox[{"Polarization", "[", 
          RowBox[{"k2", ",", 
           RowBox[{"-", "I"}], ",", 
           RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], ".", 
        RowBox[{"GS", "[", 
         RowBox[{"k2", "-", "p2"}], "]"}], ".", 
        RowBox[{"GS", "[", 
         RowBox[{"Polarization", "[", 
          RowBox[{"k1", ",", 
           RowBox[{"-", "I"}], ",", 
           RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}], "]"}], ".", 
        RowBox[{"Spinor", "[", 
         RowBox[{
          RowBox[{"Momentum", "[", "p1", "]"}], ",", "0", ",", "1"}], "]"}]}],
        "*", 
       RowBox[{
        RowBox[{"SMP", "[", "\"\<e\>\"", "]"}], "^", "2"}]}], ")"}], "/", 
     "t"}]}], ")"}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        SuperscriptBox["\<\"e\"\>", "2"], " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           FormBox["\<\"\[CurlyPhi]\"\>",
            TraditionalForm], 
           FormBox["\<\"(\"\>",
            TraditionalForm], 
           FormBox[
            OverscriptBox[
             FormBox["p1",
              TraditionalForm], "_"],
            TraditionalForm], 
           FormBox["\<\")\"\>",
            TraditionalForm]}], ")"}], ".", 
         RowBox[{"(", 
          FormBox[
           RowBox[{
            OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
            FormBox[
             RowBox[{
              OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
              FormBox["k1",
               TraditionalForm], ")"}],
             TraditionalForm]}],
           TraditionalForm], ")"}], ".", 
         RowBox[{"(", 
          FormBox[
           RowBox[{
            OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
            FormBox[
             FormBox[
              RowBox[{
               OverscriptBox[
                FormBox["k2",
                 TraditionalForm], "_"], "-", 
               OverscriptBox[
                FormBox["p2",
                 TraditionalForm], "_"]}],
              TraditionalForm],
             TraditionalForm], ")"}],
           TraditionalForm], ")"}], ".", 
         RowBox[{"(", 
          FormBox[
           RowBox[{
            OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
            FormBox[
             RowBox[{
              OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
              FormBox["k2",
               TraditionalForm], ")"}],
             TraditionalForm]}],
           TraditionalForm], ")"}], ".", 
         RowBox[{"(", 
          RowBox[{
           FormBox["\<\"\[CurlyPhi]\"\>",
            TraditionalForm], 
           FormBox["\<\"(\"\>",
            TraditionalForm], 
           FormBox[
            RowBox[{"-", 
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"]}],
            TraditionalForm], 
           FormBox["\<\")\"\>",
            TraditionalForm]}], ")"}]}]}], "t"]}], "-", 
     FractionBox[
      RowBox[{
       SuperscriptBox["\<\"e\"\>", "2"], " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          FormBox["\<\"\[CurlyPhi]\"\>",
           TraditionalForm], 
          FormBox["\<\"(\"\>",
           TraditionalForm], 
          FormBox[
           OverscriptBox[
            FormBox["p1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\")\"\>",
           TraditionalForm]}], ")"}], ".", 
        RowBox[{"(", 
         FormBox[
          RowBox[{
           OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
           FormBox[
            RowBox[{
             OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
             FormBox["k2",
              TraditionalForm], ")"}],
            TraditionalForm]}],
          TraditionalForm], ")"}], ".", 
        RowBox[{"(", 
         FormBox[
          RowBox[{
           OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
           FormBox[
            FormBox[
             RowBox[{
              OverscriptBox[
               FormBox["k1",
                TraditionalForm], "_"], "-", 
              OverscriptBox[
               FormBox["p2",
                TraditionalForm], "_"]}],
             TraditionalForm],
            TraditionalForm], ")"}],
          TraditionalForm], ")"}], ".", 
        RowBox[{"(", 
         FormBox[
          RowBox[{
           OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
           FormBox[
            RowBox[{
             OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
             FormBox["k1",
              TraditionalForm], ")"}],
            TraditionalForm]}],
          TraditionalForm], ")"}], ".", 
        RowBox[{"(", 
         RowBox[{
          FormBox["\<\"\[CurlyPhi]\"\>",
           TraditionalForm], 
          FormBox["\<\"(\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{"-", 
            OverscriptBox[
             FormBox["p2",
              TraditionalForm], "_"]}],
           TraditionalForm], 
          FormBox["\<\")\"\>",
           TraditionalForm]}], ")"}]}]}], "u"]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        SuperscriptBox["\<\"e\"\>", "2"], " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           FormBox["\<\"\[CurlyPhi]\"\>",
            TraditionalForm], 
           FormBox["\<\"(\"\>",
            TraditionalForm], 
           FormBox[
            RowBox[{"-", 
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"]}],
            TraditionalForm], 
           FormBox["\<\")\"\>",
            TraditionalForm]}], ")"}], ".", 
         RowBox[{"(", 
          FormBox[
           RowBox[{
            OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
            FormBox[
             RowBox[{
              SuperscriptBox[
               OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
              FormBox["k2",
               TraditionalForm], ")"}],
             TraditionalForm]}],
           TraditionalForm], ")"}], ".", 
         RowBox[{"(", 
          FormBox[
           RowBox[{
            OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
            FormBox[
             FormBox[
              RowBox[{
               OverscriptBox[
                FormBox["k2",
                 TraditionalForm], "_"], "-", 
               OverscriptBox[
                FormBox["p2",
                 TraditionalForm], "_"]}],
              TraditionalForm],
             TraditionalForm], ")"}],
           TraditionalForm], ")"}], ".", 
         RowBox[{"(", 
          FormBox[
           RowBox[{
            OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
            FormBox[
             RowBox[{
              SuperscriptBox[
               OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
              FormBox["k1",
               TraditionalForm], ")"}],
             TraditionalForm]}],
           TraditionalForm], ")"}], ".", 
         RowBox[{"(", 
          RowBox[{
           FormBox["\<\"\[CurlyPhi]\"\>",
            TraditionalForm], 
           FormBox["\<\"(\"\>",
            TraditionalForm], 
           FormBox[
            OverscriptBox[
             FormBox["p1",
              TraditionalForm], "_"],
            TraditionalForm], 
           FormBox["\<\")\"\>",
            TraditionalForm]}], ")"}]}]}], "t"]}], "-", 
     FractionBox[
      RowBox[{
       SuperscriptBox["\<\"e\"\>", "2"], " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          FormBox["\<\"\[CurlyPhi]\"\>",
           TraditionalForm], 
          FormBox["\<\"(\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{"-", 
            OverscriptBox[
             FormBox["p2",
              TraditionalForm], "_"]}],
           TraditionalForm], 
          FormBox["\<\")\"\>",
           TraditionalForm]}], ")"}], ".", 
        RowBox[{"(", 
         FormBox[
          RowBox[{
           OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
           FormBox[
            RowBox[{
             SuperscriptBox[
              OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
             FormBox["k1",
              TraditionalForm], ")"}],
            TraditionalForm]}],
          TraditionalForm], ")"}], ".", 
        RowBox[{"(", 
         FormBox[
          RowBox[{
           OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
           FormBox[
            FormBox[
             RowBox[{
              OverscriptBox[
               FormBox["k1",
                TraditionalForm], "_"], "-", 
              OverscriptBox[
               FormBox["p2",
                TraditionalForm], "_"]}],
             TraditionalForm],
            TraditionalForm], ")"}],
          TraditionalForm], ")"}], ".", 
        RowBox[{"(", 
         FormBox[
          RowBox[{
           OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", 
           FormBox[
            RowBox[{
             SuperscriptBox[
              OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
             FormBox["k2",
              TraditionalForm], ")"}],
            TraditionalForm]}],
          TraditionalForm], ")"}], ".", 
        RowBox[{"(", 
         RowBox[{
          FormBox["\<\"\[CurlyPhi]\"\>",
           TraditionalForm], 
          FormBox["\<\"(\"\>",
           TraditionalForm], 
          FormBox[
           OverscriptBox[
            FormBox["p1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\")\"\>",
           TraditionalForm]}], ")"}]}]}], "u"]}], ")"}]}], 
  TraditionalForm]], "Output",
 ImageSize->{467, 211},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[2]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"%", "//", 
   RowBox[{
    RowBox[{"DoPolarizationSums", "[", 
     RowBox[{"#", ",", "k1", ",", "0"}], "]"}], "&"}]}], "//", 
  RowBox[{
   RowBox[{"DoPolarizationSums", "[", 
    RowBox[{"#", ",", "k2", ",", "0"}], "]"}], "&"}]}]], "Input",
 CellLabel->"In[3]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{
    FractionBox["1", 
     RowBox[{"t", " ", "u"}]], 
    RowBox[{
     SuperscriptBox["\<\"e\"\>", "4"], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}]}]}], "+", 
   RowBox[{
    FractionBox["1", 
     RowBox[{"t", " ", "u"}]], 
    RowBox[{
     SuperscriptBox["\<\"e\"\>", "4"], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}]}]}], "+", 
   RowBox[{
    FractionBox["1", 
     SuperscriptBox["u", "2"]], 
    RowBox[{
     SuperscriptBox["\<\"e\"\>", "4"], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}]}]}], "+", 
   RowBox[{
    FractionBox["1", 
     SuperscriptBox["t", "2"]], 
    RowBox[{
     SuperscriptBox["\<\"e\"\>", "4"], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}], " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         RowBox[{"-", 
          OverscriptBox[
           FormBox["p2",
            TraditionalForm], "_"]}],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$98", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        OverscriptBox["\[Gamma]", "_"], "\[CenterDot]", "(", 
        FormBox[
         FormBox[
          RowBox[{
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"], "-", 
           OverscriptBox[
            FormBox["p2",
             TraditionalForm], "_"]}],
          TraditionalForm],
         TraditionalForm], ")"}], ")"}], ".", 
      SuperscriptBox[
       OverscriptBox["\[Gamma]", "_"], 
       FormBox[
        FormBox[
         RowBox[{"$MU", "\[InvisibleApplication]", 
          RowBox[{"(", "$96", ")"}]}],
         TraditionalForm],
        TraditionalForm]], ".", 
      RowBox[{"(", 
       RowBox[{
        FormBox["\<\"\[CurlyPhi]\"\>",
         TraditionalForm], 
        FormBox["\<\"(\"\>",
         TraditionalForm], 
        FormBox[
         OverscriptBox[
          FormBox["p1",
           TraditionalForm], "_"],
         TraditionalForm], 
        FormBox["\<\")\"\>",
         TraditionalForm]}], ")"}]}]}]}]}], TraditionalForm]], "Output",
 ImageSize->{434, 266},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[3]="]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[Cell["\t", "ExampleDelimiter"],
  $Line = 0; Null]], "ExampleDelimiter"],

Cell[TextData[{
 "This is a small piece of the matrix element squared for ",
 Cell[BoxData[
  RowBox[{
   RowBox[{"g", " ", "g"}], " ", "\[Rule]", " ", 
   RowBox[{"Q", " ", Cell[TextData[Cell[BoxData[
     FormBox[
      OverscriptBox["Q", "-"], TraditionalForm]], "InlineFormula"]], 
     "InlineFormula"]}]}]], "InlineFormula"],
 ". The proper summation over the polarizations of the gluons requires a \
choice of two auxiliary vectors (unless we subtract the unphysical \
contributions using ghosts). It is customary to take the 4-momentum of \
another gluon as the auxiliary vector in the summation formula."
}], "Notes"],

Cell[TextData[{
 "The option ",
 Cell[BoxData["ExtraFactor"], "InlineFormula"],
 " is used to average over the polarizations of the initial gluons."
}], "Notes"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ClearAll", "[", 
   RowBox[{"s", ",", "t", ",", "u"}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"FCClearScalarProducts", "[", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "p1", "]"}], "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SP", "[", "p2", "]"}], "=", "0"}], ";"}]}], "Input",
 CellLabel->"In[1]:="],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  FractionBox["1", 
   RowBox[{
    SuperscriptBox["s", "2"], " ", "SUNN", " ", 
    RowBox[{"(", 
     RowBox[{"1", "-", 
      SuperscriptBox["SUNN", "2"]}], ")"}], " ", 
    SuperscriptBox["u", "2"]}]], "2", " ", 
  SuperscriptBox[
   RowBox[{"SMP", "[", "\"\<g_s\>\"", "]"}], "4"], " ", 
  RowBox[{"SP", "[", 
   RowBox[{"k1", ",", 
    RowBox[{"Polarization", "[", 
     RowBox[{"p2", ",", 
      RowBox[{"-", "\[ImaginaryI]"}], ",", 
      RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
  RowBox[{"SP", "[", 
   RowBox[{"k1", ",", 
    RowBox[{"Polarization", "[", 
     RowBox[{"p2", ",", "\[ImaginaryI]", ",", 
      RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"2", " ", 
     SuperscriptBox["s", "2"], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k1", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}], 
    "+", 
    RowBox[{"2", " ", "s", " ", 
     SuperscriptBox["SUNN", "2"], " ", "t", " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k1", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}], 
    "+", 
    RowBox[{"s", " ", 
     SuperscriptBox["SUNN", "2"], " ", "u", " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k1", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}], 
    "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["s", "2"], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k1", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}], 
    "+", 
    RowBox[{"2", " ", "s", " ", 
     SuperscriptBox["SUNN", "2"], " ", "t", " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k1", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}], 
    "+", 
    RowBox[{"s", " ", 
     SuperscriptBox["SUNN", "2"], " ", "u", " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k1", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}], 
    "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["SUNN", "2"], " ", 
     SuperscriptBox["u", "2"], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", 
         RowBox[{"-", "\[ImaginaryI]"}], ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}], " ", 
     RowBox[{"SP", "[", 
      RowBox[{"k2", ",", 
       RowBox[{"Polarization", "[", 
        RowBox[{"p1", ",", "\[ImaginaryI]", ",", 
         RowBox[{"Transversality", "\[Rule]", "True"}]}], "]"}]}], "]"}]}]}], 
   ")"}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 FormBox[
  RowBox[{
   FractionBox["1", 
    RowBox[{"N", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       SuperscriptBox["N", "2"]}], ")"}], " ", 
     SuperscriptBox["s", "2"], " ", 
     SuperscriptBox["u", "2"]}]], 
   RowBox[{"2", " ", 
    SubsuperscriptBox["g", "s", "4"], " ", 
    RowBox[{"(", 
     FormBox[
      RowBox[{
       FormBox[
        OverscriptBox[
         FormBox["k1",
          TraditionalForm], "_"],
        TraditionalForm], 
       FormBox["\<\"\[CenterDot]\"\>",
        TraditionalForm], 
       FormBox[
        RowBox[{
         SuperscriptBox[
          OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
         FormBox["p2",
          TraditionalForm], ")"}],
        TraditionalForm]}],
      TraditionalForm], ")"}], " ", 
    RowBox[{"(", 
     FormBox[
      RowBox[{
       FormBox[
        OverscriptBox[
         FormBox["k1",
          TraditionalForm], "_"],
        TraditionalForm], 
       FormBox["\<\"\[CenterDot]\"\>",
        TraditionalForm], 
       FormBox[
        RowBox[{
         OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
         FormBox["p2",
          TraditionalForm], ")"}],
        TraditionalForm]}],
      TraditionalForm], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"2", " ", 
       SuperscriptBox["N", "2"], " ", "s", " ", "t", " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}], "+", 
      RowBox[{"2", " ", 
       SuperscriptBox["N", "2"], " ", "s", " ", "t", " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}], "+", 
      RowBox[{
       SuperscriptBox["N", "2"], " ", "s", " ", "u", " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}], "+", 
      RowBox[{
       SuperscriptBox["N", "2"], " ", "s", " ", "u", " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}], "+", 
      RowBox[{"2", " ", 
       SuperscriptBox["s", "2"], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}], "+", 
      RowBox[{"2", " ", 
       SuperscriptBox["s", "2"], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k1",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}], "+", 
      RowBox[{"2", " ", 
       SuperscriptBox["N", "2"], " ", 
       SuperscriptBox["u", "2"], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            SuperscriptBox[
             OverscriptBox["\[CurlyEpsilon]", "_"], "*"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}], " ", 
       RowBox[{"(", 
        FormBox[
         RowBox[{
          FormBox[
           OverscriptBox[
            FormBox["k2",
             TraditionalForm], "_"],
           TraditionalForm], 
          FormBox["\<\"\[CenterDot]\"\>",
           TraditionalForm], 
          FormBox[
           RowBox[{
            OverscriptBox["\[CurlyEpsilon]", "_"], "(", 
            FormBox["p1",
             TraditionalForm], ")"}],
           TraditionalForm]}],
         TraditionalForm], ")"}]}]}], ")"}]}]}], TraditionalForm]], "Output",
 ImageSize->{515, 142},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[2]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"%", "//", 
    RowBox[{
     RowBox[{"DoPolarizationSums", "[", 
      RowBox[{"#", ",", "p1", ",", "p2", ",", 
       RowBox[{"ExtraFactor", " ", "->", " ", 
        RowBox[{"1", "/", "2"}]}]}], "]"}], "&"}]}], "//", 
   RowBox[{
    RowBox[{"DoPolarizationSums", "[", 
     RowBox[{"#", ",", "p2", ",", "p1", ",", 
      RowBox[{"ExtraFactor", " ", "->", " ", 
       RowBox[{"1", "/", "2"}]}]}], "]"}], "&"}]}], "//", 
  "Simplify"}]], "Input",
 CellLabel->"In[3]:="],

Cell[BoxData[
 FormBox[
  RowBox[{"-", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SubsuperscriptBox["g", "s", "4"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}], "-", 
         RowBox[{
          SuperscriptBox[
           FormBox[
            OverscriptBox[
             FormBox["k1",
              TraditionalForm], "_"],
            TraditionalForm], "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"s", " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["N", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "t"}], "+", "u"}], ")"}]}], "+", 
            RowBox[{"2", " ", "s"}]}], ")"}]}], "+", 
         RowBox[{"s", " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["N", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "t"}], "+", "u"}], ")"}]}], "+", 
            RowBox[{"2", " ", "s"}]}], ")"}]}], "-", 
         RowBox[{"2", " ", 
          SuperscriptBox["N", "2"], " ", "s", " ", "t", " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}], "-", 
         RowBox[{
          SuperscriptBox["N", "2"], " ", "s", " ", "u", " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}], "-", 
         RowBox[{"2", " ", 
          SuperscriptBox["s", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}], "+", 
         RowBox[{"2", " ", 
          SuperscriptBox["N", "2"], " ", 
          SuperscriptBox["u", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["k2",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}], "-", 
         RowBox[{
          SuperscriptBox["N", "2"], " ", 
          SuperscriptBox["u", "2"], " ", 
          SuperscriptBox[
           FormBox[
            OverscriptBox[
             FormBox["k2",
              TraditionalForm], "_"],
            TraditionalForm], "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            FormBox[
             OverscriptBox[
              FormBox["p1",
               TraditionalForm], "_"],
             TraditionalForm], 
            FormBox["\<\"\[CenterDot]\"\>",
             TraditionalForm], 
            FormBox[
             OverscriptBox[
              FormBox["p2",
               TraditionalForm], "_"],
             TraditionalForm]}], ")"}]}]}], ")"}]}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"N", " ", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["N", "2"], "-", "1"}], ")"}], " ", 
       SuperscriptBox["s", "2"], " ", 
       SuperscriptBox["u", "2"], " ", 
       SuperscriptBox[
        RowBox[{
         FormBox["\<\"(\"\>",
          TraditionalForm], 
         FormBox[
          RowBox[{
           FormBox[
            OverscriptBox[
             FormBox["p1",
              TraditionalForm], "_"],
            TraditionalForm], 
           FormBox["\<\"\[CenterDot]\"\>",
            TraditionalForm], 
           FormBox[
            OverscriptBox[
             FormBox["p2",
              TraditionalForm], "_"],
            TraditionalForm]}],
          TraditionalForm], 
         FormBox["\<\")\"\>",
          TraditionalForm]}], "2"]}], ")"}]}], ")"}]}], 
  TraditionalForm]], "Output",
 ImageSize->{545, 141},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}},
 CellLabel->"Out[3]="]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell["", "SectionFooterSpacer"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 Cell[BoxData[
  InterpretationBox[
   StyleBox[
    GraphicsBox[{},
     BaselinePosition->Baseline,
     ImageSize->{6, 0}],
    CacheGraphics->False],
   Spacer[6]]]],
 "See Also"
}], "SeeAlsoSection",
 WholeCellGroupOpener->True,
 CellID->1255426704],

Cell[TextData[{
 StyleBox[ButtonBox["Polariazation",
  BaseStyle->"Link",
  ButtonData->"paclet:FeynCalc/ref/Polariazation",
  ButtonNote->"Polariazation"],
  FontFamily->"Verdana"],
 ", ",
 StyleBox[ButtonBox["Uncontract",
  BaseStyle->"Link",
  ButtonData->"paclet:FeynCalc/ref/Uncontract",
  ButtonNote->"Uncontract"],
  FontFamily->"Verdana"],
 "."
}], "SeeAlso"]
}, Open  ]],

Cell[" ", "FooterCell"]
},
Saveable->False,
ScreenStyleEnvironment->"Working",
WindowSize->{725, 750},
WindowMargins->{{-10, Automatic}, {Automatic, -8}},
WindowTitle->"DoPolarizationSums",
TaggingRules->{
 "ModificationHighlight" -> False, "ColorType" -> "SymbolColor", "LinkTrails" -> 
  GridBox[{{
      RowBox[{
        ButtonBox[
        "FeynCalc", ButtonData -> "paclet:FeynCalc/guide/FeynCalc", 
         BaseStyle -> {"Link", "DockedLinkTrail"}]}]}}, ColumnAlignments -> 
    Left], "HasOptions" -> False, "ExampleCounter" -> 1, "NeedPlatMsgIn" -> 
  None, "RootCaptions" -> "", 
  "Metadata" -> {
   "built" -> "{2020, 6, 20, 19, 4, 18.296011}", 
    "history" -> {"9.0", "", "", ""}, "context" -> "FeynCalc`", 
    "keywords" -> {}, "specialkeywords" -> {}, 
    "tutorialcollectionlinks" -> {}, "index" -> True, "label" -> 
    "Feyn Calc Symbol", "language" -> "en", "paclet" -> "FeynCalc", "status" -> 
    "None", "summary" -> 
    "DoPolarizationSums[exp, k, n] sums over physical (transverse) \
polarizations of an external massless vector boson with momentum k, where n \
is an auxiliary 4-vector from the gauge-dependent polarization sum formula. \
DoPolarizationSums[exp, k, 0] replaces the polarization sum of an external \
massless vector boson with momentum k by -g^\\[Mu]\\[Nu]. This corresponds to \
the summation over all four polarizations, including the unphysical ones. \
DoPolarizationSums[exp, k] sums over the three polarizations of an external \
massive vector boson with momentum k and mass k^2.", "synonyms" -> {}, 
    "tabletags" -> {}, "title" -> "DoPolarizationSums", "titlemodifier" -> "",
     "windowtitle" -> "DoPolarizationSums", "type" -> "Symbol", "uri" -> 
    "FeynCalc/ref/DoPolarizationSums"}, "SearchTextTranslated" -> ""},
CellContext->"Global`",
FrontEndVersion->"10.4 for Linux x86 (64-bit) (April 11, 2016)",
StyleDefinitions->Notebook[{
   Cell[
    StyleData[
    StyleDefinitions -> FrontEnd`FileName[{"Wolfram"}, "Reference.nb"]]], 
   Cell[
    StyleData["Input"], CellContext -> "Global`"], 
   Cell[
    StyleData["Output"], CellContext -> "Global`"]}, Visible -> False, 
  FrontEndVersion -> "10.4 for Linux x86 (64-bit) (April 11, 2016)", 
  StyleDefinitions -> "Default.nb"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "PrimaryExamplesSection"->{
  Cell[7255, 218, 386, 15, 31, "PrimaryExamplesSection",
   CellTags->"PrimaryExamplesSection",
   CellID->95617754]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"PrimaryExamplesSection", 75435, 2537}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[590, 21, 3150, 76, 53, "AnchorBarGrid",
 CellID->1],
Cell[3743, 99, 295, 11, 45, "ObjectNameGrid"],
Cell[CellGroupData[{
Cell[4063, 114, 1627, 44, 245, "Usage",
 CellID->2010081510],
Cell[CellGroupData[{
Cell[5715, 162, 738, 24, 31, "NotesSection",
 CellGroupingRules->{"SectionGrouping", 50},
 CellID->515142955],
Cell[6456, 188, 31, 0, 14, "SectionHeaderSpacer"],
Cell[CellGroupData[{
Cell[6512, 192, 101, 2, 27, "Input"],
Cell[6616, 196, 578, 15, 56, "Output"]
}, Open  ]]
}, Closed]]
}, Open  ]],
Cell[CellGroupData[{
Cell[7255, 218, 386, 15, 31, "PrimaryExamplesSection",
 CellTags->"PrimaryExamplesSection",
 CellID->95617754],
Cell[CellGroupData[{
Cell[7666, 237, 196, 6, 26, "ExampleSection",
 CellID->1253010737],
Cell[7865, 245, 147, 3, 49, "Notes"],
Cell[8015, 250, 192, 5, 45, "Input"],
Cell[CellGroupData[{
Cell[8232, 259, 531, 15, 45, "Input"],
Cell[8766, 276, 644, 25, 37, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9447, 306, 127, 3, 27, "Input"],
Cell[9577, 311, 3267, 137, 64, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12881, 453, 105, 2, 9, "ExampleDelimiter"],
Cell[12989, 457, 197, 4, 49, "Notes"],
Cell[13189, 463, 192, 5, 45, "Input"],
Cell[CellGroupData[{
Cell[13406, 472, 531, 15, 45, "Input"],
Cell[13940, 489, 644, 25, 37, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14621, 519, 127, 3, 27, "Input"],
Cell[14751, 524, 418, 17, 37, "Output"]
}, Open  ]],
Cell[15184, 544, 448, 7, 99, "Notes"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15669, 556, 105, 2, 9, "ExampleDelimiter"],
Cell[15777, 560, 121, 3, 32, "Notes"],
Cell[15901, 565, 216, 6, 45, "Input"],
Cell[CellGroupData[{
Cell[16142, 575, 531, 15, 45, "Input"],
Cell[16676, 592, 644, 25, 37, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17357, 622, 117, 3, 27, "Input"],
Cell[17477, 627, 947, 40, 60, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[18473, 673, 105, 2, 9, "ExampleDelimiter"],
Cell[18581, 677, 272, 9, 32, "Notes"],
Cell[18856, 688, 424, 12, 79, "Input"],
Cell[CellGroupData[{
Cell[19305, 704, 4226, 122, 283, "Input"],
Cell[23534, 828, 9058, 293, 232, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32629, 1126, 300, 9, 27, "Input"],
Cell[32932, 1137, 14021, 490, 287, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[47002, 1633, 105, 2, 9, "ExampleDelimiter"],
Cell[47110, 1637, 626, 13, 89, "Notes"],
Cell[47739, 1652, 161, 4, 32, "Notes"],
Cell[47903, 1658, 424, 12, 79, "Input"],
Cell[CellGroupData[{
Cell[48352, 1674, 4616, 123, 345, "Input"],
Cell[52971, 1799, 9176, 312, 163, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62184, 2116, 517, 15, 45, "Input"],
Cell[62704, 2133, 9601, 312, 162, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[72344, 2450, 31, 0, 29, "SectionFooterSpacer"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72412, 2455, 270, 12, 31, "SeeAlsoSection",
 CellID->1255426704],
Cell[72685, 2469, 367, 13, 56, "SeeAlso"]
}, Open  ]],
Cell[73067, 2485, 23, 0, 41, "FooterCell"]
}
]
*)

(* End of internal cache information *)

